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A detailed theoretical and numerical analysis of the localization length in alternating metamaterial- 
birefringent random layered stacks, under uncorrelated thickness-disorder, has been performed. Similar struc¬ 
tures have recently been reported to suppress the Brewster delocalization for p-polarized light, when "standard” 
isotropic layers (with positive index of refraction) are considered instead of metamaterial layers, providing a 
generic means to produce polarization-insensitive, broadband reflections. However, this enhancement of local¬ 
ization is valid for short wavelengths A compared to the mean layer thickness ao- At higher wavelengths, we 
recover the Brewster anomalies for p-polarized states impeding a remarkable localization of light. To achieve a 
better localization for a wider range of wavelengths, we replaced the conventional isotropic layers by negative- 
index metamaterials presenting low losses and constant index of refraction over the near-infrared range. As a 
result, our numerical calculations exhibit a linear dependence of the localization length with A (in the region 
5 < A/oo < 60) reducing the Brewster anomalies in more than two orders of magnitude with respect to the 
standard isotropic scheme at oblique incidence. This enhancement of localization is practically independent of 
the thickness disorder kind and is also held under weak refractive-index disorder. 


I. INTRODUCTION 

Broadband reflections from a disordered dielectric medium 
are a physical manifestation of the localization of light 0. It 
has been reported the design of such high-performance broad¬ 
band mirrors with a broader reflection band than periodic sys¬ 
tems with the same refractive indices dl. If a stack sys¬ 
tem could be devised to effectively localize both s- and p- 
polarized light over all angles of incidence, it would provide 
a means to obtain polarization-insensitive, broadband reflec¬ 
tions with potential applications that include waveguides and 
thermoelectric devices, among others 0. However, a fun¬ 
damental problem arises due to the so-called Brewster delo¬ 
calization 00, that is, the impossibility of achieving a to¬ 
tal localization of light at certain oblique incidences for p- 
polarization. The main reason is that the depth of penetration 
for p states increases several orders of magnitude as compared 
with the depth of penetration of s states 0 ]. This is an im¬ 
portant issue to take into account for potential polarization- 
insensitive broadband mirrors. 

Recently, Jordan et al. 1 reported a suppression of 
Brewster delocalization in stacks of alternating isotropic - 
birefringent random layers under uncorrelated thickness dis¬ 
order. Stacks containing a mixture of positive uniaxial and 
negative uniaxial birefringent layers have been claimed as ef¬ 
fective media to inhibit the Brewster delocalization for over 
all angles of incidence. This work was motivated by a previ¬ 
ous analysis of the nonpolarizing reflections from birefringent 
guanine and isotropic cytoplasm multilayer structure found in 
some species of silvery fish d. These authors mimic the 
guanine-cytoplasm structure in fish skin where there are two 
different types of birefringent crystal present with isotropic 
cytoplasm gaps and obtain a noticeable enhancement of local¬ 
ization. Nevertheless, as we will show briefly, this enhance¬ 


ment of localization for p-po lari zed light is valid for short 
wavelengths compared to the mean layer thickness. At higher 
wavelengths, we recover the Brewster anomalies for p states 
at oblique incidence, impeding a remarkable localization of 
light. 

In order to achieve an improvement of localization, 
we propose a stack composed of alternating metamaterial- 
birefringent disordered layers where the isotropic layers pos¬ 
sess a negative index of refraction (that is, we are dealing 
with left-handed (LH) isotropic media ifTTl fl2h ). These ar¬ 
tificial materials have been reported to resolve images be¬ 
yond the diffraction limit lfl3lll4ll . act as an electromagnetic 
cloak [HI, enhance the quantum interference |[r7.] or yield 
to slow light propagation [JL8]. Moreover, recent research 
on LH metamaterials formed by hybrid metal-semiconductor 
nanowires |19}] has described negative indices of refraction 
in the near-infrared with values of the real part well below 
-1 and extremely low losses (an order of magnitude better 
than present optical negative-index metamaterials). With this 
choice, we demonstrate a reduction of Brewster anomalies 
in more than two orders of magnitude with respect to the 
standard isotropic scheme at higher wavelengths, providing 
a means to obtain polarization-insensitive, broadband reflec¬ 
tions in a wider wavelength region. Furthermore, this en¬ 
hancement of localization is practically independent of the 
thickness disorder kind and is also held under weak refractive- 
index disorder. 

The plan of the work is as follows. In Sec. II we carry 
out an exhaustive description of our system and the theoret¬ 
ical framework to derive an analytical expression for the lo¬ 
calization length £ SiP in the short-wavelength regime (that is, 
when the incident wavelength A is less than the mean layer 
thickness ao). A detailed numerical analysis concerning the 
localization length for two different disorder distributions (the 
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Figure 1. Our disordered layered stack of alternating uniaxial posi¬ 
tive (n+) and uniaxial negative (n_) layers. Gray regions correspond 
to left-handed isotropic layers with negative index of refraction no. 
The principal axes of the birefringent layers are aligned to x and y 
directions, while light enters from the left at angle 9 and propagates 
through z axis. 


Poisson-thickness disorder and the uniform-thickness distri¬ 
bution) will be performed in Sec. III. Two wavelength re¬ 
gions will be analyzed in detail: the previously mentioned 
short-wavelength regime and an intermediate-wavelength re¬ 
gion (5 < A/ao < 60) where, as we will show, the improve¬ 
ment of localization is more pronounced. For completeness, 
weak refractive-index disorder (with random variations in the 
indices of refraction of each layer) is also considered in our 
model. A set of analytical expressions for £ s p in the so-called 
intermediate-wavelength region, as a function of the wave¬ 
length A and the incident angle 9, are also derived. Finally, 
we summarize our results in Sec. IV. 


II. SYSTEM DESCRIPTION AND THEORETICAL 
MODEL 


Our structure of interest is formed by a mixture of alternat¬ 
ing uniaxial birefringent layers and LH isotropic media (see 
Fig. □>• For positive uniaxial layers, the corresponding re¬ 
fractive index vector is given by n + = (m, n\, 712 ) while for 
negative uniaxial = ( 712 , 712 , 711 ). The principal axes of 
birefringent layers are aligned to x and y directions and the ra¬ 
tio of negative uniaxial to the total number of layers has been 
defined via the parameter /. Each birefringent layer is sand¬ 
wiched between two left-handed isotropic layers with negative 
index of refraction no and the whole systems is embedded in a 
’’standard” isotropic media with positive refractive index \no |. 
The following condition 712 > rii > |no | is hold to avoid crit¬ 
ical angles, that is, to be sure that light is able to access all 
angles of incidence at each isotropic-birefringent interface in 
the stack. Light enters from the left at angle 9 and propagates 
through z axis. 

In order to effectively improve the localization of light in 
such mixed structures, we need LH metamaterials with ex¬ 


tremely low losses and constant values of the negative-index 
over a relatively wide spectrum. Such metamaterials have 
been recently reported by Paniagua-Dominguez et al. iH 
and might operate in the near-infrared region (from 700 to 
1900 nm in this case). As we will show in the next section, the 
enhancement of localization in our system is mainly achieved 
in a region where the incident wavelength A satisfies the in¬ 
equality 20 < A/ao < 60, where ao corresponds to the mean 
layer thickness. Accordingly, our mixed disordered stack can 
adequately be tailored to operate in the near-infrared provided 
that ao — 35 nm. 

Let us now perform a theoretical study of the localization 
length for both s and p states in such structures under uncor¬ 
related thickness disorder. To this aim, we derive an expres¬ 
sion for the transmission coefficient of the whole arrangement 
r r s ]J and then we use the standard definition of the localization 
length 

t _ 2 L 

S S,p — — Tj Tj (1) 

(In T S 'p) 

where L = Nao is the system length and N the total number 
of layers. The angular brackets (...) stand for averaging over 
the disorder. The theoretical method here employed is a non- 
perturbative approach based on the exact calculation of the 
Green’s function (GL) of an electromagnetic wave in a given 
electric permittivity profile 120). Once the GL of the layered 
structure is known, the transmission coefficient T s p can be 
derived using the Lisher-Lee relation [21]. This GF approach 
is compatible with the transfer matrix method and has been 
widely used to calculate the average density of states over a 
sample, the energy spectrum of elementary excitations [22] or 
the characteristic barrier tunneling time (231], among others. 

In our Green’s function approach, the inverse of the trans¬ 
mission coefficient through our multilayer structure, 1/T S)P , 
coincides with the diagonal elements of the resulting transfer 
matrix and can be derived via the following relation [20] 

T s ,p = \D n \~ 2 , (2) 

where Dm corresponds to the so-called characteristic deter¬ 
minant 


N 


D n = exp [~i<p N ] ]^[(1 - 1/2 . 


( 3 ) 


z=i 


Here ipM represents the total phase thickness accumulated by 
the incident wave along the whole sample, and is given by 


N 


VN 




( 4 ) 


1=1 


where tpi is the phase thickness of each individual layer. For 
our birefringent media, ipi satisfy @] 


f 2jL 
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) aini'X^Jl- ( n °^ n 0 J, p-polarization 
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with ai the layer thickness and nij (j = x, y, z) the corre¬ 
sponding indices of refraction of the positive or negative uni¬ 
axial layers 


ni,x — — 


rii, positive layers 
n 2 , negative layers, 


and 


ni, z = 


712, positive layers 
?ii, negative layers. 


( 6 ) 


(7) 


In our LH isotropic media, the parameters <pi are calculated 
via the following expression 


= 


Cbl \ tiq | COS0, 


( 8 ) 


for both s- and p-polarized states. Notice that the phase thick¬ 
nesses have negative sign in this case, due to the negativeness 
of the indices of refraction. 

I) 1 !,- 1 corresponds to the determinant of a tridiagonal matrix 
and satisfies the following recurrence relation 

D\ 0) = A,Dl°\ - B,D (0) 


'l — 

with the initial conditions 


i-ii 


(9) 


Ax = 1, D £ 0) = 1, and D ( °[ = 0. 


For l > 1 we have 

Ai = 1 


n-i ,i 

H-2J-1 


exp [2 iipi], 


and 


Bi = (A t - 1)(1 - rf_ x ;_ 2 ), 


( 10 ) 


( 11 ) 


where the quantity r;_^ symbolizes the Fresnel reflection 
amplitude of light propagating from region l — 1 into region 
l. In our case, these coefficients are given by the generalized 

m 


isotropic-birefringent Fresnel relations 

|no| cos 9 — yjnfy — ? 1 q sin 2 9 


n~i,i = 


n ly - n l sin 2 6 


|?1o| COS0 - 

for s polarization and 

\no\ni, x ni tZ cos 0 - nlJnf z - nl sin' 


( 12 ) 


n-i,i = 


\no\ni, x ni, z cos 9 + n^Jnf z - Uq sin“ 0 


(13) 


If the Fresnel coefficients rf l _ 1 are small, one can easily 
find (after simple inspection of the recurrence relation Eq. (0) 
that the parameter D^ = 1. Therefore, the determinant I)\; 
(and, consequently, the transmission coefficient r f s , p ) may be 
calculated at arbitrary distribution of the layers thicknesses. 
Provided that our system presents thickness disorder only (so 
the reflection coefficients rp_ 1 ; remain invariable along the 
averaging process) and for a large number of layers N, we 
can write 


N 

n 

i=i 


(l-rf-i ,)- 1/a 


(1 -r*) N+ (l~r> ' N ~ 


1 - 1/2 


S— ,P— z 


(15) 

where N± stands for the total number of positive or negative 
uniaxial layers, and the reflection coefficients r s± iP± are given 
by Eqs. (IT2l) and (IT3l) . with the appropriate choice of uniaxial 
layer (that is, positive or negative) and incident polarization. 
After some trivial algebra, we obtain the following analytical 
expression for the localization length £ s p (see again Eq. i|7} ) 


2 do 

£s,p 


= - In 




(16) 


where / = N-/N. Notice that, under the low-reflection as¬ 
sumption, our Eq. (fl6] > is independent of the thickness disor¬ 
der kind and the incident wavelength. We will analyze later 
the validity of this theoretical result for short wavelengths 
compared to the mean layer thickness, that is, A < do- 

After simple inspection of Eq. (IT6l >. one observes that the 
Brewster anomalies for p-po lari zed states only occur when 
/ = 0 (isotropic-positive stacks) or / = 1 (isotropic-negative 
stacks). In such particular cases, the right-hand side of Eq. 
m becomes zero at resulting in the logarithmic diver¬ 
gence of the localization length and, consequently, in the 
delocalization of p states. This result can be easily evidenced 
in Fig. [2] where we represent the inverse of the localization 
length for p-polarization, t;- 1 (calculated via Eli. ( I I (if ) versus 
the parameter / and the angle of incidence 9. Two differ¬ 
ent values of the positive birefringence A + = n -2 — n-| have 
been selected (with our choice of the refractive index vec¬ 
tors, the following relation holds A = —A + ). The index 
of refraction of the left-handed layers has been chosen to be 
no = —1.33. One observes that for real mixtures of positive 
and negative uniaxial layers |§] (that is, when 0 < / < 1) no 
Brewster anomalies occur. For higher birefringence parame¬ 
ters (see Fig. 0b)) we obtain a smoothly varying landscape for 
Consequently, as we increase the strength of the layers 
birefringence, the more localized states we achieve. 


III. NUMERICAL RESULTS 


for p-polarized states. As it is well-known, these reflec¬ 
tion coefficients for p polarization vanish at the generalized 
anisotropic Brewster angle 9b I25il 


tan(0e) 



(14) 


In this section we present our numerical results concern¬ 
ing the localization length for two different disorder distribu¬ 
tions, that is, the Poisson-thickness disorder and the uniform¬ 
thickness distribution (the latter also under weak refractive- 
index disorder). Two wavelength regions will be analyzed in 
detail: the short-wavelength regime (A/ao < 1), where we 
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Figure 2. Inverse of the localization length (calculated via Eq. 
jl6H versus the parameter / and the angle of incidence 9. Two dif¬ 
ferent values of the positive layers birefringence A+ = ri 2 —ni have 
been chosen. 



will check the validity of our main theoretical result Eq. ( IT6] i 
and an intermediate-wavelength region (5 < A/ao < 60) with 
an enhancement of localization up to two orders of magnitude 
with respect to standard isotropic layers. 


Figure 3. Localization length versus the angle of incidence 9 in 
the short-wavelength regime (A/ao = 0.032) for a Poisson-thickness 
distribution. Top figures stand for s-polarized states while bottom 
figures represent p-polarized light. Solid lines symbolize our numer¬ 
ical calculations whereas dashed lines show our theoretical results 
(Eq. GD). 


A. Poisson-thickness disorder 


Let us first consider the well-known Poisson thickness dis¬ 
tribution for our disordered system, where the averaging pro¬ 
cess of In T s p is carried out via the following expression 


1 f L , 

a 

— da exp 


«o Jo 

a 0_ 


InT, 


s,p> 


(17) 


where, for each single realization, the numerical computation 
of the quantity 7’ s p is obtained via the exact recurrence re¬ 
lation Eq. To this aim, we show in Fig. [ 7 ] the local¬ 
ization length ^s.p versus the angle of incidence 9 in disor¬ 
dered stacks with mixing ratios of / = 0.25 (left column) 
and / = 0.75 (right column) in the short-wavelength regime 
(A/ao = 0.032). Dashed lines correspond to our theoretical 
results based on Eq. ( IT6l ). One observes a good agreement be¬ 
tween our numerical and theoretical results with an error of 
less than 5%. As can be noticed from Fig. [3] in the short- 
wavelength regime the localization length of p states remains 
constant in a wide range of incident angles (about 45 degrees 
for / = 0.25 and 70 degrees for / = 0.75). As previously 
commented, higher values of the positive birefringence A + 
yields to an enhancement of localization. This is a manifesta¬ 
tion of a total suppression of Brewster delocalization, in good 
agreement with recent numerical results jjj. 

For higher wavelengths the situation is markedly different, 
as can be easily observed in Fig. [4] Here we represent our nu¬ 
merical calculations for the localization length £ s ,p as a func¬ 
tion of the wavelength A for two different angles of incidence, 
9 = 0 and 9 = 28 degrees. Let us first concentrate on the top 
figures |4ja) and (4j b ) where isotropic layers with index of re¬ 
fraction |n 0 | = 1-33 have been considered. The birefringence 
of positive layers has been chosen to be A + =2.33 and the 
ratio / = 0.75. Solid lines symbolize our numerical results 


performed via Eq. O. One notices that £ s ,p has a quadratic 
dependence on A at higher wavelengths, a typical result for 
various disordered models and the subject of intensive ongo¬ 
ing research ||2(||27]. Dotted lines correspond to the following 
numerically-deduced equations for the localization length as 
a function of the wavelength 


6 

a 0 


1-306 

a 0 


+ 


1 


1.15C 2 r 2 (cos9) 0 - 5 


(18) 


§p 

a 0 


1.306 


A 


ao 1.15C 2 r 2 (cos 9) 2 - 0 \a 0 


P 0 


(19) 


where C = 1.15(ni + ri 2 + |no|) and £ So (£po) represents the 
localization length of s ( p ) states in the limit of short wave¬ 
lengths (see again Eq. (IT6l >). In addition, r 2 p stands for the 
average interfacial reflection coefficient in our system nm 
and is given by 


r 


2 

s,p 


(! - f) r 


2 

s+,P+ 


+ fr 


s — iP— ' 


( 20 ) 


An in-deep numerical analysis has been performed to derive 
Eqs. m and ( | 1 9k To this aim, 20 different layered stacks 
at different angles of incidence have been considered, and a 
detailed study of the coefficient characterizing the quadratic 
dependence on A has been carried out. Notice that the results 
derived via Eqs. ( IT8| > and ( | 1 9| i are in good agreement with 
our numerical calculations for angles of incidence up to 28 
degrees. At higher angles of incidence (not shown in our fig¬ 
ures) the deviation between numerical and analytical results is 
greater than 5%. 

Let us now replace the ’’standard” isotropic layers by left- 
handed isotropic metamaterials (this change, as previously 


yunits of a 0 ) yunits of a 0 ) 
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mentioned, assumes that the index of refraction no is now neg¬ 
ative). As it was recently shown by our group for isotropic lay¬ 
ered structures & this change of sign results in a significant 
modification of the transmission coefficient T s p . As a conse¬ 
quence, the localization length £ SiP presents a region of linear 
dependence with A when the positions of the layer boundaries 
are randomly shifted with respect to ordered periodic values. 
The same behavior is found for LH isotropic-birefringent lay¬ 
ered stacks. 

The physical explanation for this effect is related to the van¬ 
ishing of the total phase accumulated along the structure. This 
vanishing only occurs for stacks of alternating right-handed 
media (isotropic or birefringent) and left-handed isotropic lay¬ 
ers, where the alternating phases possess opposite sign. Pro¬ 
vided strong uncorrelated thickness disorder, the total phase 
accumulated over the whole structure is practically null. Thus, 
the transmission coefficient T s p (which strongly depends on 
the total phase along the structure) has a smoother dependence 
on the frequency w than the typical right-handed stacks, where 
the average of lnT SiP (see the denominator of Eq. £[])) gets the 
well-known dependence W 2 . In our case, this average has a 
linear dependence with the frequency (that is, the localization 
length is proportional to the wavelength A). 

Our numerical calculations confirm this feature as can be 
noticed in Figs. [4}c) and [4}d). Solid lines symbolize our nu¬ 
merical calculations (Poisson thickness distribution, Eq. ( ITTb ) 
while dotted lines correspond to our numerically-deduced set 
of equations 


6 _ U , 2 

CjrsKcosfi 1 ) 0 - 5 

2 / A_ 

ao a 0 ' C|r p |(cos0) 30 \ao 

After direct comparison of Figs. HJc) and [4jd ) (left-handed 
isotropic-birefringent stacks) with Figs.[4ja) and[4jb) (’’stan¬ 
dard” isotropic-birefringent stacks) one observes an enhance¬ 
ment of the localization for isotropic metamaterial structures 
at higher wavelengths, mainly due to a smoother dependence 
of the localization length £ S}P with the incident wavelength A. 



a o a o 





B. Uniform-thickness disorder 

Let us now study another type of thickness disorder where 
the arrangement of layers has random boundaries [:2Q], For 
completeness, weak refractive-index disorder has also been 
included in our model, where we have randomly changed the 
indices of refraction of each layer via the following relations 
no = rig 0 ) +<Jj5 (isotropic layers) and n± = n±) +<Jj5 (bire¬ 
fringent layers, either positive or negative). The superscripts 
(0) indicate unperturbed refractive indices whereas the param¬ 
eters <jj are zero-mean independent random numbers within 
the interval [—0.5, 0.5], The strength of the disorder is mea¬ 
sured by the parameter S which has been chosen to be roughly 



Figure 4. Localization length £ S , P versus the incident wavelength 
A for two different angles of incidence, 9 = 0 and 6 = 28 de¬ 
grees. Top figures represent ’’standard” isotropic-birefringent stacks 
(with positive indices of refraction) while bottom figures stand for 
LH isotropic-birefringent disordered structures. The birefringence 
of the positive layers has been chosen to be A+ = 2.33 and the ratio 
/ = 0.75 in all cases. Solid lines represent our numerical calcu¬ 
lations (Poisson thickness distribution, Eq. G3) whereas the dotted 
lines symbolize the localization length calculations performed via 
Eqs. ( 1 181 and ( 1191 ) (top figures) and Eqs. OH and ( 1221) (bottom fig¬ 
ures). 


a 2% of the corresponding unperturbed indices of refraction 
(that is, we assume weak refractive-index disorder). 

Firstly, the short-wavelength regime will be considered. In 
Fig. 0we show the localization length £ s p versus the angle 
of incidence 6 in such disordered stacks with mixing ratios 
of / = 0.75 (left column) and / = 0.25 (right column) 
for A/ao = 0.032. Dashed lines correspond to our theo¬ 
retical result, Eq. (IT6t . Even though this equation was de¬ 
rived for thickness-disorder only, a good agreement between 
our numerical and theoretical results is found, provided weak 
refractive-index disorder. 

As we increase the incident wavelength, the behavior of 
the localization length is significantly different (as previ¬ 
ously stated in the Poisson-thickness distribution). In Fig. 
[6] we represent three-dimensional plots of the localization 
length £ SiP for a combined scheme of uniform-thickness 
and refractive-index disorder. Left-hand figures display the 
numerical calculations for isotropic-birefringent stacks with 
| no | = 1.33, whereas right-column graphs stand for our left- 
handed isotropic stacks with no = —1.33. In order to per¬ 
form a reliable average of the logarithm of the transmission 
coefficient, lnT s p , 800 random configurations of our disor¬ 
dered stack have been considered in all cases. Once more, the 
birefringence of positive layers was A .. = 2.33 and the ratio 
/ = 0.75. It can be noticed that, at long wavelengths, a de¬ 
crease of one order of magnitude for £ p at normal incidence 
is achieved (compare Fig. [6jc) to[6jd)) and even up to two or¬ 
ders of magnitude at higher angles. Similar numerical results 
have been obtained taking into account uniform-thickness dis¬ 
order alone (not shown in our figures). Consequently, un- 


^units of a 0 ) 
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Figure 5. Localization length £ S:P versus the angle of incidence 6 
in the short-wavelength regime (A/ao = 0.032) for an uniform¬ 
thickness distribution. Left figure stands for s polarization while 
the right one represents p-polarized light. Solid lines symbolize our 
numerical calculations whereas dashed lines show our theoretical re¬ 
sults (Eq. GD). 
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Figure 6. Enhancement of localization, under a mixed scheme 
of uniform-thickness and weak refractive-index disorder, for LH 
isotropic-birefringent stacks. Numerical calculations for £ 3:P versus 
the wavelength A and the angle of incidence 9 are shown. The bire¬ 
fringence of positive layers was chosen to be A+ = 2.33 while the 
ratio / = 0.75 in all cases. Left-column figures represent our nu¬ 
merical results for isotropic-birefringent stacks with |no| = 1.33 
whereas right-column graphs stand for our left-handed isotropic 
stacks with no = —1.33. The simulations were averaged over 800 
random configurations of the layered stacks. 


der a strong uncorrelated positional disorder scheme, weak 


refractive-index disorder does not affect significatively the lo¬ 
calization of light. 

IV. DISCUSSION AND CONCLUSIONS 

A detailed theoretical and numerical analysis of the local¬ 
ization length, under uncorrelated thickness-disorder, in alter¬ 
nating metamaterial-birefringent random layered stacks has 
been performed. Our work was motivated by the recently re¬ 
ported suppression of Brewster delocalization in alternating 
isotropic-birefringent stacks d. where Jordan et al. mimic 
the guanine-cytoplasm structure in fish skin formed by two 
different types of birefringent crystal present with isotropic 
cytoplasm gaps §]. Nevertheless, this enhancement of lo¬ 
calization for /./-polarized states is achieved for short wave¬ 
lengths (compared to the mean layer thickness do). At higher 
wavelengths, we recover the Brewster anomalies for p states 
at oblique incidence, impeding a remarkable localization of 
light. 

Our disordered metamaterial-birefringent stacks can ef¬ 
fectively improve the localization of light in the so-called 
intermediate-wavelength region (in our case, 5 < A/ao < 60) 
up to two orders of magnitude with respect to the standard 
isotropic scheme at oblique incidence. This is mainly due to 
the smoother dependence of the localization length with the 
incident wavelength (linear dependence) in this particular re¬ 
gion, unlike the well-known quadratic dependence in positive- 
index disordered stacks. Moreover, this enhancement of lo¬ 
calization is practically independent of the thickness disorder 
kind and is also held under weak refractive-index disorder. 

In order to reduce losses in our system (principally due 
to the presence of left-handed layers) we considered meta¬ 
materials with extremely low losses and constant values of 
the negative-index over the near-infrared region (from 700 to 
1900 nm in this case) in. As shown in Sec. Ill, the predom¬ 
inant enhancement of localization in our system is achieved 
in the region 20 < A/ao < 60. Accordingly, our mixed 
disordered stack can effectively be tailored to operate in the 
near-infrared provided that o.q ~ 35 nm. 

Our group consider that these metamaterial-birefringent 
disordered structures might provide great insight for the ex¬ 
perimental design of polarization-insensitive, broadband re¬ 
flectors. 
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